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Abstract. An outstanding problem in Earth science is understanding the 
method of transport of magma in the Earth's mantle. Two proposed methods 
for this transport are percolation through porous rock and flow up conduits. 
Under reasonable assumptions and simplifications, both means of transport can 
be described by a class of degenerate nonlinear dispersive partial differential 
equations of the form: 

<h + {<t> n )z - (<T«>~">t)z) z = o, 

where <f>(z,0) > and t) — > 1 as z — > ±oo. Although we treat arbitrary n 
and m, the exponents are physically expected to be between 2 and 5 and and 
1, respectively. 

In the case of percolation, the magma moves via the buoyant ascent of a 
less dense phase, treated as a fluid, through a denser, porous phase, treated as 
a matrix. In contrast to classical porous media problems where the matrix is 
fixed and the fluid is compressible, here the matrix is deformable, with a viscous 
constitutive relation, and the fluid is incompressible. Moreover, the matrix is 
modeled as a second, immiscible, compressible fluid to mimic the process of 
dilation of the pores. Flow via a conduit is modeled as a viscously deformable 
pipe of magma, fed from below. 

Analog and numerical experiments suggest that these equations behave akin 
to KdV and BBM; initial conditions evolve into a collection of solitary waves and 
dispersive radiation. As cf> — > 0, the equations become degenerate. A general local 
well-posedness theory is given for a physical class of data (roughly H 1 ) via fixed 
point methods. The strategy requires positive lower bounds on rf>(z, t). The key to 
global existence is the persistence of these bounds for all time. Furthermore, we 
construct a Lyapunov energy functional, which is locally convex about the uniform 
porosity state, = 1, and prove (global in time) nonlinear dynamic stability of 
the uniform state for any m and n. For data which are large perturbations of the 
uniform state, we prove global in time well-posedness for restricted ranges of m 
and n. This includes, for example, the case n = 4, m = 0, where an appropriate 
uniform in time lower global bound on <j> can be proved using the conservation 
laws. We compare the dynamics to that of other problems and discuss open 
questions concerning a larger range of exponents, for which we conjecture global 
existence. 
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Figure 1. (a) A single Olivine crystal, and the channels along which 
molten rock flows, 1 . (b) Fluid flow up a viscously deformable pipe. 

1. Introduction and Overview 

An outstanding problem in Earth science is understanding the method of transport 
of magma in the Earth's mantle. One proposed method is buoyant ascent of the melt 
through a viscously deformable porous medium. Such a model may be visualized in 
Figure ^ ( a ) > where the magma flows through tubules along the grain boundaries of 
the solid rock. Equations for this process were derived independently in and EJ ■ 
In both derivations, the system is modeled using a volume averaged two phase flow, 
one for the melt and one for the matrix, at an extremely low Reynolds number. This 
amounts to having two coupled conservation of mass equations to two force balance 
equations, 

d t {pf4>)+V- {pf<pv s ) = 

at(p.(i-0))+v-(p.(i-0K) = o 

V • (<f>a { ) = Pf<pg + I 

V-{{i-<f>)a s ) = Ps {i-4>)g-l 

The variables are defined in tabled The fluid rheology is that of an incompressible 

Table 1. Viscously Deformable Media Variables 



Parameter 


Symbol 


Porosity 




Fluid density 


Pf 


Solid density 


Ps 


Fluid velocity 


Vf 


Solid velocity 


V s 


Fluid rheology 


°~f 


Solid rheology 


OS 


Gravitational vector 


9 


Interphase force 


I 



inviscid fluid while the solid rheology is treated as viscously compressible fluid. We 
refer the reader to [3 El El El E] for expositions on these relationships and their 
simplifications. 

When reduced to 1-D, as in El El III El E| j this transport can be described by 
a class of degenerate nonlinear dispersive partial differential equations of the form: 



{<p m - dA n B z ) \r m M = -n (i.i) 

This can be rewritten as the coupled system, 

4>t - ct> m u (1.2) 

(r i -d z rd z )u = -(r) z (1.3) 
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The exponent n corresponds to the power in the permeabily-porosity relationship 

K cx (j> n 

while the exponent m corresponds to the power in the bulk viscosity-porosity 
relationship 



1 




In 0], the authors conclude that the parameter space for these exponents is 
2 < n < 5 and < m < 1. In [2], the bulk viscosity is taken to be constant, 
corresponding to m = 0, and 2 < n < 3. Several papers [HI 003 have focused on the 
case n = 3 and m — 0. 

A second important mechanism for magma migration is that of flow through a 
viscously deformable pipe embedded in a viscous matrix, as takes place in the thermal 
plumes of the convecting mantle. This is pictured in Figure 0(b). The primitive 
equations, [TT)| . 

d t A = - d z Q 

Q = (Apg + mi d z {A~ l d z Q)) 

whose variables are defined in Table |2 give rise to an equation of the form , 
with n = 2 and m = 1, replacing <f> with A. Analog laboratory models for 

Table 2. Conduit Flow Variables 



Parameter 


Symbol 


Areal region of pipe 


A 


Flux up pipe 


Q 


Fluid viscosity 


VL 


Matrix viscosity 




Fluid-Matrix density difference 


Ap 


Gravitational constant 


9 



this mechanism were studied in |1(JI 111), verifying the theory by demonstrating the 
appearance and interaction solitary waves, with measurements consistent with the 
equations. Throughout this paper, we will use <j) as the dependent variable. 

Upon inspection, it is not clear what the behavior of 1)1.1)1 is. It has both 
a nonlinear forcing term and a nonlinear dispersive term, in addition to being 
degenerate and nonlocal. The degeneracy is evident in 1)1.3)) . where the inevitability 
and elliptiticity of the operator break down as <fi — > 0. This mathematical challenge is 
at the heart modeling the underlying physical system; can the melt disaggregate in the 
case of porous flow, or can the conduits split in the case of pipe flow? Are there upper 
and lower bounds for physical variables such as porosity? 

Previous studies, [S1[7|EJE1^SE1E1E1i 01 the equations did not successfully 
address the issue of well-posedness and the related issue of a lower bound on <j>, the 
(scaled) melt fraction, although ^5] did point out several reasons why solutions for 
which cf> went to zero would be non-physical. 

Degenerate nonlinear parabolic partial differential equations are an important 
class of equations arising in porous media flow (with a fixed matrix) ^B] and in 
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geometric flows, e.g. motion of a surface by its local curvature, |19j . However, 
the study of degenerate nonlinear dispersive equations appears to be wide open. 
Our equations bear some resemblance to those considered in [20], known to admit 
compactly supported solitary waves. We can see the dispersive character of the 
equation by linearizing about some background porosity <f>Q. The linearized equation 

4> t - ncj)^ 1 - (f>Q~ m <l) Z zt = 
has the dispersion relation 



n<f>Q k 
1 + ^~ m k 2 



u{k) -- 
with group velocity 

W ' (fc) " (l + <j> n - m k 2 ) 2 
When (f>o — 1, which will be the background value of <f> in our analysis, we see the 
similarity between this equation and the Regularized Long Wave (RLW) equation, 
also known as the Benjamin-Bona-Mahoney (BBM) equation 

d t u + d x u + ud x u — d^dtu = 
which has a linearized (about zero) dispersion relation 

" (fc) = TTV 2 

with group velocity 

1 - k 2 

<"'(*) = (iTW 

The BBM equation |^ |22 , was proposed as an alternative to the Korteweg - de 
Vries (KdV) equation and is asymptotically equivalent to KdV in the regime of small 
amplitude long waves. Like KdV, BBM admits stable solitary waves, although it is 
not completely integrable. The magma equations also have solitary waves, which play 
a central role in the general dynamics; see, for example, pfl IH31 151 111 l*Hfl 171 IT21 ITU 151 19] . 
In this paper we focus on well-posedness questions and stability of the uniform state. 
The question of stability of solitary waves is currently under study. 
The main results of this paper are the following: 

(1) A local well-posendess theory is established for arbitrary initial conditions, <p(z, 0), 
in a physically natural set of functions, X , of finite "energy" , £. The time-interval 
of existence, and therefore global well-posedness, is shown to be controlled by 
|| 0|| h 1 + min ze R (f>(z, t). See Theorem 12. 121 

(2) We obtain ranges of material exponents, (m, n), for which global well-posedness 
for arbitrary data in X holds. See Corollaries 14.61 and 15.71 along with figure |2 
This range of parameters includes a set of obviously degenerate cases (such as 
m = and n = 4), as noted in Corollary 15.81 and discussed in Section |SJ The 
important case of McKenzie's equation, m = and n = 3 remains open. We 
conjecture global well-posedness. 

(3) For any values of m and n, we prove that the uniform porosity state, (/>(z, t) = 1 is 
shown to be nonlinearly stable in Corollary 14. 71 The key step is to show that an 
appropriately constructed functional, defined on X, is locally convex about the 
uniform state. This strategy has been used in many Hamiltonian and nonlinear 
wave problems, see |2H I2ll2l)l HE\ . 
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Figure 2. The trapezoid shows the permeabily-viscosity pairs of exponents, 
(n,m), for which we have global existence for data of arbitrarily large H 1 norm. 
The line n + 3m = 4 bounds the region for which we can ensure a lower bound, 
while the line n + 3m = 6 bounds the region for which we can ensure an upper 
bound. The line n + m = 3 is a technical constraint. For all other (n, m) pairs, 
we have global existence for data near the uniform state of 1. 



In a forthcoming paper, we shall treat the question of stability of small amplitude 
solitary waves [77] , 

We proceed as follows. In section [2 we exhibit local well-posedness in H k (M.) 
spaces for k — 1,2,... In section [3] we discuss certain conserved quantities associated 
with the equations. Section |3| addresses certain properties of a conserved functional 
and the nonlinear stability of the uniform background state. In section 0] we prove 
global well-poesdness for data near the uniform background state, and in section [5] 
we show global well-posedness for large data for a certain subset of the equations. In 
section El we make comparisons with other equations, conjecturing how the equation 
regularizes itself through the time dependent length scale of the problem. We then 
summarize and discuss future work and open problems. 

Throughout this paper, we shall use the following notation. 

(i) H k denotes the Sobolev space W k ' 2 (R) of functions with L 2 weak derivatives up 
to order k. k will always denote an integer, k > 1. 

(ii) C k,a denotes the space functions with derivatives up to order k are Holder 
continuous with exponent a 

(iii) All integrals are taken over R in the spatial variable or on some finite interval 
[0, t] in the temporal variable. 

(iv) Generic constants are denoted by C, C , C" ... or D, D',. . . 

(v) The dependence of a solution <f> of 1(1.1(1 on its data be may be explicitly expressed 
as 

0(z,i;0 o ) 
where (f>o is t ne initial condition. 

2. Local in time well-posedness in H k , k > 1 

In this section we consider the well-posedness of the initial value problem for (|1.1|) . 
existence, uniqueness, and continuous dependence with respect to the data. Formally, 
solving 1(1. 3|) for it, we substitute this solution into 11.2(1 and integrate with respect to 
t. Applying an appropriate initial condition for 0(z,t), we obtain 

m = ^o + [ 4> m L^[-{<t> n )z]ds, (2.1) 
Jo 

L^u = - (<b n u z ) z + cj) m u (2.2) 

The mapping A is highly nonlinear, in no small part due to the dependence of 
the nonlocal operator, L^, on <f>. Our strategy is to construct solutions of the initial 
value problem for 1(1.211.3(1 seeking a fixed point of the mapping <fi i— > A[</>] in an 
appropriate metric space of functions. In particular, we shall apply the contraction 
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mapping principle 28 to A, restricted to an appropriately chosen closed subset of 
functions, Xj^ e T , for which 4>(z,t) — 1 G H 1 . The set Xj^ eT is constructed so that 
iteration of A on it preserves ellipticity and invertibility properties of La. 



Definition 2.1 (i) Given 1 > e > 0, R,T > 0, fc > 1, we define 

l R,e,T 



,*)-l||ir* <^> supH^.,*)- 1 ]^ <e- 1 }(2.3) 



sup 

t<T t<T 

The sets X^ € T obviously include the constant solution, <fi — 1, so they are 
nonempty and closed sets. 

(it) A metric on X^ eT is defined by 

[</> - ^] k {T) = Sup U(; t) - 1>(;t)\\ Hl , (2.4) 
t<T 

We note that A" is a metric space, not a function space; it is centered about the 
constant solution 1. 

Definition 2.2 We call 4> a local solution of with data 4>o on the time interval 
[0,T], if(f>>Q,(f>-leH k and(j> = A[<p}. 

Theorem 2.3 (Local well-posedness) Let 1 > e > 0, R > 0, and k > 1 . Then, for 
T > sufficiently small, A is a contraction on X^ T with respect to the metric !T\ ) 
and hence admits a unique fixed point. Thus, the initial value problem for fJ.gl 7^ ) 
has a unique local solution. Moreover the solution depends continuously on the initial 
data. 

Remarks on notation: Generic constants, C, may depend on R,e,T and k, along 
with the two parameters in the equation n and m. While the dependence on n and 
to will suppressed in our proofs, the reader should be aware of it, and the change in 
the behavior of C as e \ and R /* oo as n and to vary. In the cases of geophysical 
interest, where n > 2 and < to < 1, constants may behave like: 

C — > oo, as e \ and R /* oo 

Thus, ensuring that cj> stays away from zero is critical to our problem as it evolves in 
time. In the following, unless otherwise specified, k > 1, R > and 1 > e > 
are arbitrary. We now embark on setting up the proof of Theorem l2.3l which requires 
several technical steps. 

Proposition 2.4 Let R, e, T > and k > 1. For all <f> £ X^ £ T 

(i) 4> e C* -1,1 / 2 (R) 

(ii) (f> > e 

(Hi) There exist Mj < oo such that < Mj for j < k — 1. The Mj depend on 

R, e and k. 

Proof: This follows from Sobolev inequalities [2U and the definition of X^ e T . 

□ 

Also of use is that the Sobolev spaces H k (R) form an algebra and admit the following 
inequality give citation 
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Lemma 2.5 Assume k > 1/2. If f,g G H k (R), then 

\\fg\\ Hk <C k \\f\\ Hk \\g\\ Hk (2.5) 

The nonlinearity in the magma equations appears as a power of the dependent variable. 
It will be useful to note the following property of these functions. 

Proposition 2.6 Let 1 > e > 0, R,T > and k > 1. For all p £ K there exists a 
Lipschitz constant C = C(R, e, k,p), such that for 0, ip € X^ T , 

IW-.iip-VManiH* <C||<K-,ii)-V>(-,*2)|| H * 

/orh.ia < T. 

Proof: This follows from the bounds of Proposition l2~H Lemma |2~B1 and that property 
that x h- > cc p is C°° for a; bounded away from zero. 

□ 

Proposition 2.7 (i) Let R,e,T > and fc > 1. Let <E e T and f £ L 2 . Then, 
L<t>(-,t)U = f has a unique solution u £ iJ 1 (R) for t £ [0, T}. Moreover, there 
exists a constant C — C{R, e), such that 

IK-,*)|| ffl <C||/|| i2 (2.6) 

(ii) Furthermore, if cp £ X^ + e 2 T and f £ H k , then u(-,t)H k+2 £, then there exists a 
constant C = C(R, e, fc) such that 

\\u\\ Hk+2 <C\\f\\ Hk (2.7) 
Proof: See j2H] and for details. We outline the arguments. 

(i) Since <f> £ T , L^ is self-adjoint and positive definite. The bilinear form 
(f/,V% = f4> n UVdz + / (j) m d z Ud z V dz defines inner product on H 1 . By 
the Riesz Representation Theorem and there exists a unique H 1 solution, u(-, t), 
which satisfies (V, u(-, = JVf dz for all V £ H 1 . Clearly (£5) holds. 

(ii) Higher regularity may be deduced by first proving it for the case fc = and 
then applying induction. The fc = case is proved by noting that the regularity 
of (j> implies the coefficients will be C k+1 and studying the limit of difference 
quotients to estimate higher derivatives. Alternatively, since the problem is in 
one spatial dimension, we could compute and estimate derivatives of u explicitly 
using variation of parameters representation of the L^ 1 . 

□ 

Proposition 2.8 Let 1 > e > 0, R, T > and k > 1 > 1. There exists a constant 
C = C(R, e, fc) such that for all tj> £ tT , if u solves L^u = —(4> n ) Z) then 

\H H *<c 

Proof: The result follows by applying Proposition 12 . 71 to / = ~(<fi n ) z £ H k ~ 1 . 

□ 

Proposition 2.9 (i) Let l>e>0, R> and fc > 1. Assume \\<fro — l\\ Hk < R and 
0o > 2e. Then there exists T\ — T\(R, e, fc, (j>o) > 0, such that for T <T% , 

<f>£*L,T — A[<t>]£X k (2.8) 
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(ii) Moreover, there exists < Ta < T 1; Ta = T 2 (i?, e, A;, </>o) > such that for T < T 2 
A is a contraction on T , i.e., there exists an a — a(R, e, k, <f>o) < 1 such that 
for all <f),ip € Afjj c T , 

[A[0]-A[0]] t (T)<a[0-^(T) (2.9) 

We first prove Proposition ^. 9l (|i)l. We must thus find a Ti such that choosing T < T\, 
for all * < T 



UA^G,*)-!!!^ <i? 



A[«/>](,i) 



< 



and 



A[^](.,*)-i>ic 1 ([o J r 1 );fr*) 

To establish (|2~TU|) . let u = ^I 1 [-(</>")*]. Then 



(2.10) 
(2.11) 

(2.12) 



-11 



\ci) m (;s)u(;s)\\ Hh ds 



Applying Lemma \'2. 51 Proposition 12 . 61 and Proposition 12 . 71 to the integrand we get 
W n u\\ Hk <CU m - l\\ Hk \\u\\ Hh + \\u\\ Hk < C"{CC M-lWv +1)<C (2.13) 
Hence, 

||A0-1|| H , < \\1> -I\\ m +CT 

Choosing || <po — l\\ Hk < R, T may be chosen sufficiently small such that 1 1 Acf> — l\\ Hk < 
R. This shall be our first candidate for T%, which we denote as t\. 
To prove 1)2. 11[) . consider the difference between </>o and A<j>, 

l|A0(-,t) - ML < V^||A^(-,t) - Mm < V2 / H m (;s)u(;s)\\ H1 ds < V2CT 

Jo 

where C is the constant in (|2.13|) . 

\m\ > \m - \m - m > \m - iia[0] - ^oiioc > 2e - v^ct 

Taking T sufficiently small, this will be bounded below by e. Hence, 



i 



Af 



1 

< - 

e 



l|2.12|l is more difficult to establish, requiring the following two Lipschitz-type 
estimates, which will also be needed in the proof of Proposition ^. 91 (fTTf) . This additional 
work is necessary because of the appearance of nonlinearity in the operator itself in 
contrast to more common nonlinear nonlocal dispersive equations. We thus requires 
bounds on the "closeness" of the operators for different 0's, which we establish in 
the following propositions. 

Proposition 2.10 Assume l>e>0,T,R>0 and k > 1. Let <f>, ip be in X^ eT . 
Then for any two values t x , t 2 < T at which <f> and ip are evaluated, there exists a 
constant C — C(R, e, k) such that 
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Proof: Let u and v satisfy L^. M) u = -((/>(■, h) n ) z and L^. t t 2 )V = -(V>(-> t 2 ) n ) z . 
Then 

L<j>[u] ~ Ly[v\ = ~(cf> n u z ) z + {il> n v z ) z + <Tu - ^ m v = {i> n - <t> n )z 
Let 7 be in H (K). Multiplying by 7 and integrating by parts, 



n u z - rv z )i z + {<Tu - v m «)7 = <rh 

After adding and subtracting and rearranging terms, 

+ (V m - <£">7 (2-14) 
For fc = 1, if we substitute u ~ v for 7, we get the estimate 

ran - v) z f + r\u - vf <w- <n L2 n« - »ii La 

+ n^-0"ii^ii Vz ii i2 ii( M - U ) z ii i2 

+ ||^ n -0 n || L o O |M| ia ||«-t;|| ia 

which, after some further manipulation becomes, 

"Hhi ^ c \\<t> ~ ^Ww 
For greater regularity, with k > 2, we return to l|2.14|l . and integrate by parts 



r(u - «) 2 7, + 4> m {u - «) 7 = / - <TR), + r - <W" - <T» 7 

This expresses u> = u — v as a solution to the elliptic problem L^w — / with 

/ = -((v>" - 4> n ) Vz ) z + r- - <T> 

Indeed, for k > 2, we know w is at least a weak solution, as / 6 L 2 , and the coefficients 
in the operator are bounded. We may thus apply Proposition 12 .71 to get 

\\u-v\\ Hk <C\\f\\ Hk - 2 <CU-4,\\ Hk 

□ 

We also need to be able to consider normed differences of <j) m L^ > 1 [—((f) n ) z ] for which 
we make the following estimate, 

Proposition 2.11 Assume 1 > e > 0, T, R > and k > 1. Given (f> i/> in A^ er , 
i/iere exists a constant C = C(i?, e, A;) smc/i i/iai /or any t\, t<i < T 

Proof: Let u and v satisfy L^u = —(4> n ) z and L^v — —(i/j n ) z evaluated at t\ and t 2 , 
respectively. Then, using Propositions 12.61 |2~B1 and 12.101 

\\4, m u~^ m v\\ Hk < U m u-q> m v\\ Hk +U m v-yj m v\\ Hk 

< (c u m i\\ H h + 1) n« - v\\ Hk + c \\ v \\ Hk u m - f n \\ Hk 
<cu-n Hk 

□ 
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We are now able to prove (|2.12|) . First we will establish differentiability, and then 
continuity, of A [</>]. Let h ^ 0. Using Proposition ^. Ill 



AM(.,i + /0-Afe](.,t) 



*)«(-,*) 



H k 



^ +h <r{;s)u{;S)-^{;t)u{;t)ds 



< 



< 



\h\ 

c_ 

\h\ 



l+h 



t+h 



\^ m (; S )u(; S )-r(;t)u(;t)\\mds 



Using C 1 continuity of <j> in time, there exists r\ such that 

||<H-,s)-#,t)IU < 77 1« — *j 

thus, 

AM;t + h)-A[<P](.,t) 



So as /i — > 0, 



AM(-,t + fc)-AM(-,t) 



ft 



^(•,t)u(.,t) 



0"*(.,t)u(.,t) 



< H 



t+h 



Is — i| ds 



c4 



(2.15) 



Thus it is differentiable, hence continuous, and 

dtM4>](;t) = 4>{-,t) m L^ t) [-((/>(-,t) n ) z ] 
Now, let t 1 ,t 2 < Ti. Using Proposition |2~TD and (l2~T5l) . 

||W(-,t 2 )-cW(-,ii)|| H * = ||0" l (-,i 2 )u(-,i 2 )-0 m (-,t 1 ) w (-,t 1 )|| H1 

<C-||<K-,t 2 )-#,fi)|| ff * 

Since is continuous in t, we are done. 

We have thus established that for T < T\ , (|2.12l) , and hence Proposition 12.91 
holds, where T\ is determined from the values necessary for 12.10fl and (|2.11() to hold. 

Now we will establish that A is a contraction. Let t <T\. Applying Proposition 

l2~TTl 



Jo 



ds 



< / cu-^w^ds^ct^-^it) 

Jo 

If we now choose t such that Ct < 1, say T 2 , then we have established 
[A0-AV>]k < a[0-V]ft 

With 

a = CT 2 < 1 

We have thus established that for data, </>o satisfying conditions 
\\<j>o-l\\ ah <R 



(2.16) 



(2.17) 
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and 



< 



2e 



(2.18) 



for some i?>0,l>e>0, and k > 1. There will exist a T2 > such that the 
map A will be a contraction on the space Xj^ e Ta proving Proposition ^. 91 (uj). By the 
contraction mapping theorem, there exists a unique fixed point in X^ T , proving 
existence and uniqueness in Theorem 12. 31 Reformulated as a local existence theorem 
for (|l.lfl . we have 

Theorem 2.12 Let 4>q satisfy 



II 



< R 



< 



a 4> - 1 e 



2e 

for R>0, 1 > e > 0, and k > 1. 

There exists a T > dependent on R, e and R — 
C^flO.T) : H k (R)) such that 

<K>,t) = <h + J cj ) (-,s) m u(-,s)ds 

-(#•, t) n u{; t) z ) z + 0(;t) m u(;t) = -(#, t)"), 

T/izs is a local solution to $1.11 such that <fi(-,t) > e and ||0(-,t) — < Rfort < T. 
Moreover, there is a maximal time of existence, T max > 0, such that i/T max < oo, 

then 

1 



lim 

t-» T ma 



(2.19) 



Proof: The first part of this theorem has already been proved. Let us address the 
notion of the maximal time of existence. Assume T max < 00 and 

1 



sup 

0<t<T ma > 



< K < 00 



Then we have \\(j){-,t) — l\\ H h < K and (/>(■, t) > K^ 1 . Hence we may apply our local 
existence theory with R = K + 1 and 2e = K^ 1 to get a time of size <5 = S(K, k) over 
which we know we may propagate forward data satisfying these constraints. If we 
then take as a new initial condition (f>i(z) = 4>{z, T max — 6/2), we know it will continue 
up till at least time T max + 6/2. This contradicts the maximality of T max , proving 

i2~m 

□ 

The next result establishes that the solution to the initial value problem of l|l. lfl 
depends continuously on the data 4>q. 

Theorem 2.13 Let (j>o and ipo both satisfy 

1 



< 



2c 
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for some l>e>0,i?>0 and k > 1, and let 4> and ip be their respective solutions in 
j, for some T > that both solutions are known to satisfy. Then there exists a 
constant C = C(R, e, k) such that 

U(-,t)-^{-,t)\\ Hk < I^o-V'oIIh* e ct 

for t < T. 
Proof: Let t < T. 



U(;t)-^(;t)\\„ k < 







< \\4>o ~ ^o\\ H " + / \\M-, S ) - ^t(-,s)\\ Hk ds 
Jo 

Applying Ij2.15|l and Proposition 12 . 1 1 f withf = ti = £2), 

U(;t)-^(;t)\\ Hk < Ho-Mm+C l U{-^)-^{;s)\\ Hk 

Jo 

GronwalPs inequality then gives 



□ 



We have thus shown that the problem is locally well-posed in time for appropriate 
data, completing the proof of Theorem 12. 31 

3. Conservation Laws 

In this section we discuss conservation laws of (11.111 . We derive a particularly useful 
"energy" , which we denote by £ and then prove the continuity of the mapping <j> t— *• £ [(/>] 
on a subset of H 1 . Finally, we show that £ is locally convex about the uniform state, 
0=1. These observations play a central role in the coming discussion of global 
well-posedness and stability in sections |H and [5j 

A systematic search for conservation laws of (|l.lfl was carried out in [31] . For 
general pairs, (n, m), two conserved quantities were found. The first, corresponding 
to conservation of mass, is 



T x = J (j)-ldz (3.1) 

is conserved. The second quantity, T2 depends on the particular (n, m) pair, forming 
the following three families 

' 2rn (j) 2 z + log (j) - (j) + ldz for n + m = 1 



2 T 

A = l I i</>~ 2m <^ - log cj)dz for n + m = 2 (3.2) 

1 A2-n-m _ 1 

-r 2m ^ + r~ - Tw j_ Txdz for n + m ? 1, 2 

2 (n + m — l)(n + m — 2) 

Let T2 be the integrand of T^. Although T2 is the same within each of these three 
families, the corresponding flux, Z2, satisfying <9 t T2 + d z zi — differs depending on 
whether or not m = 1. See |3*T1 for the full taxonomy. 
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It is interesting to note that neither 7\ nor T 2 ar ^ well-defined on the spaces in 
the space, where local well-posedness is proved! In particular, (j> — 1 e H k does not 
imply Ti[4>] < oo. Nevertheless, we have found that an appropriate linear combination 
of T\ and T 2 is well-defined and can be used in the study of global well-posedness. 

This observation is linked with question of stability of the uniform state, = 1. To 
facilitate our study of the stability of the uniform state, we seek a linear combination 
of T\ and T 2 , such that: 

(a) £[4>] = a\T\[(j)\ + a 2 T 2 [</>], is a continuous functional of and well-defined on local 
solutions 

(b) <fi = 1 is a critical point of £ [0] 

(c) £[4>] is locally convex at = 1. 

We find that these criteria, (a)-(c), can be satisfied if we choose ct\ and a 2 as follows: 

for n + m = 1 

11 for n + m = 2 ,„ n 

ai = < (3.3) 




for n + m 7^ 1,2 

a 2 = 1 (3.4) 



Our functional then has the formf 

£W\ = J \r 2m 4>l + V n , m (4>) = aiTi + a 2 T 2 (3.5) 

where the potential V nym function depends on the particular pair of exponents, and is 
given as 

x log x — x + 1 for n + m = 1 

, . , , . - log x + x - 1 for n + m, = 2 

K,m(x) = < (3.6) 
,2-n-m_ l + f n + m _2)f x - I) 

— '-- n + m^l,2 

[n + m — l)(n + m — 2) 

Let us make a few remarks on V njln . 

(i) V n , m is nonnegative and C°°([e,M]), with e > and M < 00. 

(ii) K n , m (l) = V^ m (l) = 0, hence for p < 2 

(x - i)p 

is also in C°°([e,M}). 

(iii) V^i.m is continuous as a function of n and m. Indeed, if e = n + m — 2, then 

. , , x' e - 1 + e(x - 1) , x _c - 1 

lim V„ m (x) = lim ; = lim h x — 1 = — log x + x — 1 

<r^o ' e ^o (e + l)e e-»o e 

Similarly, ife = n + m— 1, 

hmK )m (x) = hm — 

x 1_e — x 

= lim (x — 1) = x logx — x + 1 

e^O e 

I In the case n = 3, m = 0, £ plays a role in the study of the instability of solitary waves in higher 
dimensions with respect to transverse perturbations 0. 
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Hence there is really a single expression not only for the potential but also for 
the functional, and the boundary cases n + m = 1 and n + m = 2, can be derived 
from it by taking limits. 

Proposition 3.1 The functional £ is Lipschitz continuous on the set 

{</>- 1 € H\R) : ||0- 1|| H1 < > e} 
under the H 1 norm, with Lipschitz C — C(R,e,n,m), i.e. 

\S[^}-S[^]\<CU-i;\\ H1 

Proof: Let <j> and tp be two elements of this space. Let M = 1 + V%R, which bounds 
by a Sobolev inequality. Since it will be of use, let I denote the Lipschitz 
constant for V n , m (x) / (x — l) 2 on the set [e, M] and let N denote this quotient's L°° 
norm on this set. Then, 



\£[<t>]-£[ 



l 

< - 

~ 2 
1 

< - 
- 2 



1 

< - 
- 2 



i> 2 



+ / \V n ,m{4>) ~ Vn,m{lp)\ 



<t>l- 


V> 2 


4> 2m 


Ki,m(0) 


(0- 


- 1) 2 


1 






/ 




{4>- 


-I) 2 



1 



1 



1p 2 



1)2 _k^a) (Vj _ 1)2 
-1/>x\\4z + 1l>z\ + 



(j) 2m - i> 



2m 



(f) 2m 1p 2 



-1) 2 -(V^-1) 2 |+ / \(Tp-l) 2 \\V n , m (<t>)-V n , m {iP)\ 

<c\\<t> z - ip z \\ L2 + c u z - ^n ffI + c" u - n L * + c" u- 

□ 

The following result shows that </>„ = 1 is a critical point of £ and for any (n, m), 
£ is locally convex near <j>„ . 

Theorem 3.2 Let p E H^R). Then 

£l±+p} = l\\p\\m+0(\\p\\m) 

as WpWw -> °- 

Proof: Each of the cases of (|3.6|) must be examined separately, however the scheme 
is universal. We present the proof for n + m > 2. Given p, let / = 1 + p. The first 
variation of £ [4>] is 



(5£[l+p],p) 



1 



ft 



n + m -I n + m-lf n+m - 1 f 1+2 

, _ ft 



1 



1 



^n + m-1 n + m - 1 /l+2m 
And the second variation is 

1 



fzz 



(S z £[l+p]p,p) 



f 2 f 
m(l + 2ro) .„,'„_ + 2m 



f 



l+2m 



P 



1 

7^ 



Jn+m ■■~\- ■ ■■-/ j2+2m 

Thus, taking a Taylor Expansion about p = 0, 

£[1 +p] = £[1] + (6£[l],p) + l(S 2 £[l]p, P ) + 0(\\pf Hl ) = ~ \\p\\ 2 Hl + 0(b|| 3 Hl ) 
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4. Data near the uniform state: Global existence and Lyapunov stability 

We now prove global well-posedness for data sufficiently close to the uniform state, 
= 1. Furthermore, we show Lyapunov stability of the uniform state. An outline of 
our strategy is 

(i) Establish conservation of £ [0] for solutions of H 2 spatial regularity. 

(ii) Construct H 1 a priori estimates using conservation of £[0] and Theorem lH . 21 when 
is of H 2 spatial regularity 

(hi) Construct H 2 a priori estimates by deriving a growth estimate on ||0 zz || i2 for 
H 2 solutions to (|l.l|l 

(iv) Use these bounds to get H 2 global existence 

(v) Exhibit stability in the H 1 norm for H 2 solutions with ||0g — l\\ H i sufficiently 
small 

(vi) Establish that £ is also conserved for H 1 solutions by approximating H 1 data by 
H 2 data and taking appropriate limits 

(vii) Prove H 1 global existence by using the H 2 argument, with the H 2 norm replaced 
with H 1 now that we have £[4>] = £[4>o] 

Proposition 4.1 Assume 0- 1 € C 1 ([0, T) : if 2 (R)) is a solution to then £[(j)} 

is conserved. 

Proof: Since is of H 2 spatial regularity, we may expand l|l.lfl to 

<t> t + {cj> n + m^^Mz ~ r- m <f>tz)z = 

Multiplying by 

' fl - 6 1 - m - n ) 



n + m — 1 

for n + m =/= 1 and log0 for n + m = 1, one will find, after integrating by parts, the 
expression may be rearranged into 

d t e + d z f = 

where, for n + m ^= 1,2, 

_ 1 02 02-m-n _ 1 + ( n + m _ 2) (0 - 1) 

6 ~ 2 2m (n + m- l)(n + m-2) 

/ = — — ^ { 0" + m0"-™- 1 t ^ - 0"" m tz - T ^0 1 " m - m0 2m z t + ^- 2m ^ 

n + m — 1 \ 1 — to 

and J e = £. Analogous expressions hold in the other cases. Note that the factor by 
which we multiply is continuous in the limit of n + to = 1 . 

□ 
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Proposition 4.2 Assume 0—1 £ C 1 ([0,T) : iJ 2 (]R)) is a solution to and there 

exist constants -R, e, smc/i i/iai 

H(;t)-1\\<R 



< 



for t < T. Then there exist constants Co = c$(R, e, n, m) and c\ = Ci(R,e,n,m) such 
that 

-||<M^)ll 2 <^ + ci||<M->*)ll 2 

Proof: Expanding out Equation <|1.1() . 

b n - l 4> z + m{n ™ i\-An-m-i^a 
Which we rewrite as 



t + np- ~<p z + m(n - m - l)^"™- 1 ^* - (n - 2m)0"- ro - 1 z zt + mf-^V^zz - 0"~"> t 



n- 



Letting u satisfy L^it = 

It 



m(n — m — 1) — - (n — 2m) — h m — - — 



l )z, <frt = </> m w, and 

4>lu 



>zzt 



- i • ' - " 4> 2 z u 4> z u z uo z 

— t; h n—, h m(n — m — 1 — ^ to ?i — 2to — ^ n — 2m — : h to—; — 



Now, let us study the evolution of the ||0zz|| hi time 
rfl 
eft 2 



d 1 IU II 2 
3Io ll^lla 



< 



'■>zz<Pzzt 

1 



Z2l - — ~ h n-n h to(to — 1)-t; (n — 2m)—; h 



An — 2m 



+ |to(1 — TO 



Il0**ll 2 Hl 2 + M 

I 



1 



||<MI 2 ll<M 



A2-m 



ll</MI 2 ||^ M || 2 + l n ~ 2m l 



||0zz|| 2 ||0z«z|| 



ll<MI 2 

< c || U || 3 1|^|| 2 + C u zz \\ 2 + c" H^IL iklL ||0«|| a 
+ c""ll^lloolKII 2 ll^ll 2 + c""Nlooll^ll 2 

Recall Proposition 12 .81 that there exists a constant depending on R, e, and k = 1 such 



that 



f/ 1 



< C 



Also, an application of Sobolev's inequality gives the bound 

WMao ^ ^Uz\\ 2 Hzz\\ 2 < V2RU ZZ \\ 2 

Using these bounds, 



^||<Ml2 < c oHzz\\ L 2 + C 1 \\(t> z z\\ 2 L2 



Therefore, 



dt 



\\</>zz\\ 2 < C + Cl ||0 

zz || 



□ 
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Figure 3. The two curves provide bounds on £<j> for \\<j> — ^\\ H i < i/2/4. In this 
picture we have assumed that the peak of the lower bound is less than %/2/4, hence 
r) c is the peak. Taking r\ < rj Ct we are then able to constrain the energy by choosing 
100 — llliji < <5 to £[<p(-, t)] < E. This assures us that ||<^(-,t) — l||jji < rj. 



Corollary 4.3 Under the assumptions of Lemma \4-<\ 

||d 2 0(-,i)|| 2 <c te Cl * + ||c^ || 2 e Clt 
Proof: Apply Gronwall's inequality. 

□ 

We note that no assumption on the size of the data was used in proving 
Propositions 14. Il and l4. 21 they will also be used in the case of "large" data in sectional 

Theorem 4.4 ( Global existence in H 2 and stabilty near the uniform state in H 1 ) 

There exists r/ c = rj c (n, m) > such that for any rj < r\ c there exists 6 = S(r], n, m) 
such that for any </>o satisfying 

\\<j> - < oo (4.1) 
Uo-M\ H t<S (4.2) 

there exists a unique <f>, <f> - 1 S (^([OjOo) : H 2 (R)) such that solves fT7|) with 
0(-,O) = 0o- Moreover, 

H(;t)-l\\ H r <V 

for all time. 

Proof: First, let us construct r\ c . Proposition 14.11 ensures that £ is conserved for 
such a solution. From Theorem 13.21 there exist constants C and D such that for 
H-l\\ Hl <V2/4, 

l\\4>-M\ 2 H1 -C\\<t>-lf H1 < £[4>] < \ U - lf m + D ||0- l\\ 3 H1 (4.3) 

The peak of l/2x 2 - Cx 3 occurs at x = 1/(3C), and we set r] c = min{l/(3C), \/2/4}. 
Let E c = l/2r] 2 - C?] 3 . Given n < n c , let E = l/2r/ 2 - Crf and let S solve 
E = 1/2S 2 + DS 3 . As indicated in figure this construction forces 

S < V < Vc < V2/4 
and for \\(f> — l\\ H i < 6 

E[4>] <E<E C 

Now assume that ||0o — l\\ H i < 5 in addition to having ||0o — M\h 2 < 00 • F rom 
Theorem EU2 we know there exists T > such that - 1 G C x ([0, T) : H 2 (R)) that 
solves (|1 . 1|> with 0o as the initial condition. 

Suppose for some time T, — 1 ceases to be in C 1 ([0, T) : iJ 2 (]R)). Furthermore, 
let T be the minimal such value. If T — oo, then we are done, so we may take T < oo. 

For t < T, - 1 € C^QM : # 2 (K)); hence it is a strong solution to (fTTTft and 
£[0] is conserved. This implies 

||0(-,t)-l|| ffl <^ 

for t < T. 
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Suppose not. Then there exists t < T for which the reverse holds, and by the 
continuity of 0, there exists a minimal time t\ such that r\ < ||0(-, < r/ c . 
However, 

/2 

U(;h)\\ H x <e c < Y 

So 1)4. 3|l holds with the same constants C and D, but this implies, 

m-,ti)} > ~ ||0(., *i) - IHhx - C||^(-,ti) - l||^ > \r! 2 - Crf = E> £[cj> Q \ 

contradicting conservation of £. 

By construction rj < r\ c < \/2/4, so (f> > \. Applying Proposition ^. 21 with R = r\ c 
and e = 5, we have that for t < T, 



; t) -l\\ H 2 <r, + c Q te c ^ + U(-,0) z z\\ e Cl * 
V2 



: ^ + c Q Te c ' T + \\<f>(-,0) zz \\e^ T = R 



< 2 



(4.4) 
(4.5) 



Hence for all t < T, 

U(;t)\\ H > 



<K-, t) 



< R + 2 < 00 



and this sum does not go to infinity, as required by l|2.19|l . for the solution to have a 
finite time of existence. Therefore, the solution persists for all time. Additionally, we 
have 



,t)-i\\ Hl <V 



□ 



Proposition 4.5 Let (j) satisfy gl| and let (f>(z,t;(/> ) - 1 € C^^T) : H^R)) be 
the solution of the initial value problem of with this data. Then £[4>{', t; 0o)] is 
conserved on the interval [0,T). 

Note that T > by the assumptions on the data and Theorem 12. 121 
Proof: Let {</>o,n} be a sequence satisfying 



II 



H 1 



< OO 



||0o,™ - < 8 
in addition to assuming that 



lira 



«0 - <P0,n\ 



H 



1 =0 



(4.6) 



Here, we have used the density of H 2 within H 1 to find such functions. 

Since </>o,« have global solutions by Theorem l4.4l we can use Theorem l2.13l on the 

space X 1 1 , to get 



\\<t>(-,t;(f>Q,m) - <A(-,*;<Ao,n)|lHi ^ ||^0,m - <^0,n|| ff i e 

Letting m — > 00, for t < 00, 



(4.7) 
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Now, consider the change in the functional £, over time, 

\£[<t>{;t; M] - £[<M\ < \£[<t>(; *5 M] - £[<K;t] <h,n)]\ 

+ \8[<l>(.,t; <h, n )] ~ £[<h,n]\ + \£[<h,n] ~ £[4>o]\ 

Since ip(z,t;ipo >n ) are global strong solutions, £[<j>(z,t;4>o,n)] is conserved. Applying 
Proposition 13 . II with the set 

W -i e ^( ffi ):||0-l|| H1 < e ,0>i} 

we get 

\£{<j>{;t; M] - £[M\ < C(U(-,t; cf> ) - 4>(-,t; 4> , n )\\ H1 + Ho - <t>0,n\\ H r) 
Now, using l|4.7|l . 

|£[0(-,t;0o)]-£fa>]| < C ]|0 O — 00,^11^ 
Letting n — > oo, we have for t < T, 

£[^{;t;^)]=£[M 

□ 

Corollary 4.6 Theorem \4-4\ holds for H . 

While we have not proved global existence in C 1 ([0,T) : _ff fc (]R)) for k > 2, the 
results of this section still yield a stability result in these smoother spaces. We know 
we have a local solution in these spaces from Theorem 12. 121 Proposition ^. II and 14.51 
establish conservation of £ for solutions in C 1 ([0,T) : H k (M.)), with k > 1. Following 
the same argument as used in theorem 14.41 we thus have 

Corollary 4.7 (Nonlinear Stability of the Uniform State) There exists e c > such 
that for any e < e c there exists S such that if \\<f>o — l\\ H i < 5, then the solution to 
fT7)) . (j> G C^flOjT) : ff fe (K)), fc > 1, T > 0, allowing T — oo, with this initial 
condition satisfies 

for t < T. 

5. Global Well-Posedness for (almost) arbitrary data 

We now turn our attention to the case of data that may not be close to the uniform 
state in the H 1 norm. While the data must still be bounded away from zero, we are 
able to show global well-posedness for data of arbitrarily large size in the H 1 norm, 
provided we satisfy the following three conditions on (n,m). 

n + m > 3 (5.1a) 
n + 3m>4 (5.16) 
n + 3rn < 6 (5.1c) 

This is the shaded region in figure The requirement that m > for smaller values of 
n is consistent with our suspicion that having a variable bulk viscosity regularizes the 
problem, as this requires progressively more work to be done to expel fluid from the 
pores as the porosity becomes smaller. For the remainder of this section, we assume 
these conditions to hold. 
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We must derive certain conserved quantities in this special case. Recalling (|3.5[) . 
we write 

g r (x) = rx r+1 - (r + l)x r + 1 (5.2) 
With r = n + m — 2, we may write 

1 4>l 1 grQ) 

2 2m r(r + 1) <\f 
Lemma 5.1 For r > 1, x > 0, 

g r {x)> {x-lf 

Proposition 5.2 Assuming £[</>] < oo and </> > 0, there exists M = M(£[<fi],n,m) < 
oo, such that 



Proof: 



< M 



-(d>{z)-lf = I {4>( Z )-l)cf> z {z)dz< I nil)- l\\n,\,lz 



\$z\ I'/'- 1 ! m+(n+m-2)/2 
Am J ) {n+m-2)/2 ™ 



< 



i(n+3m-2)/2 



1 



(ra+m-2)/2 



We are able to move the exponent outside the L°° norm since n + 3m — 2 > by 
l)5.1b(l . Then employing (|5. 1 af) to use Lemma 15. II 

< u\t +3m - 2)/2 vMr + i)m 



Hence, 



or 



Moo - 1 < 11^ - llloo < ^ ll^ll^ +3m_2)/4 (2r(r + 1)) 1/4 V^M 

M| l-(»+3m-2)/4 < ^ (2r(r + 1))V 4 y^J + ! 

which ensures an upper bound so long as 1 — (n + 3m — 2)/4 > 0, which is precisely 

□ 

This will ensure is pointwise bounded from above in terms of the data through 
conservation of £. Let us define the upper bound functional that corresponds to this 
maximum as 



u[4>] 



l-(„+3m-2)/4 



^V / 2(2r(r + l)) 1/4 v /£p] + l 



(5.3) 



Proposition 5.3 Assuming £[(j)] < oo and (j) > 0, there exists e — e{£ [</)], n,m) > 
such that 



1 



< 



1 



Degenerate Dispersive Equations 
Proof: 



21 



log- 



1 



— dz = — 



i+r/2 



l-(t> 4 



-dz 



oo 



T ) r i+r/2 ^dz 



< 



< 



1 




0, 




1-0 


-/" 




oo 


(j) m 


2 


<^l 2 


2 J — oo 



i+r/2-1 



y/2r(r + T)£ [ft + 1 

< W[0](" +3m - 4 )/V2r(r + l)£[</)] + 1 = log - 

where we have explicitly used l|5.16l) to move the exponent outside of the L°° norm, 
and we have implicitly used the other two conditions along with this one, in order to 
make use of the conserved upper bound ,U. 

□ 



We denote the corresponding functional for e as C, given as 

£[0] = exp {-W[0] (n+3m - 4)/ V2r(r + l)£[ij>] - l} 
This will ensure a conserved lower bound exists by conservation of £ . 



(5.4) 



Proposition 5.4 Assuming £[(j>\ < oo and<f> > 0, there exists R = R(£[<f>],n,m) < oo 
such that 

Proof: We divide our proof into two cases (i) : m > and (ii) : m < 0. For m > 0, 



£\4>\ > 



l 



l 



> min 



2U[(j>] 2m z U[c/)} n + m - 2 r{r + 1) 
1 1 



ifdz 



2U[<f>] 2m ' U[<j>] n + m - 2 r(r + 1) 

and we have the constant R. 
For m < 0, 

CUV 2 ™- ,n 1 



m > 



> min 



2 Yz U[4>] n+m - 2 r(r + 1) 
£\6}- 2m 1 



2 ' U[(f>} n+m - 2 r{r + 1) 



Lh 1 



l?dz 



□ 



Again, let us define a corresponding functional that provides a conserved bound on 
the H 1 norm. 



l£\ 



n<t>\ = { 



i 



i 



2U[<j)] 2m ' U[<j)] n + m '~ 2 r{r + 1) 



' £ \d>] min 



i 



for m > 



for rn < 



(5.5) 



2 ' W[0] n + m - 2 r(r + 1) 

Now that we have a priori bounds on <f> pointwise from above and below, together 
with a bound on the H 1 norm, we may prove global existence. Thus we have, 
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Theorem 5.5 Let cf> satisfy, 

Ho 
1 



1\\ H 2 < OO 



< OO 



(5.6) 
(5.7) 



Then there exists a unique <t>, <j>- 1 € ^([O.oo) : # 2 (R)) such that solves ifTTT)) 
0(-,O) = <^o- Moreover, []1/^>|| < oo /or time. 

Proof: This proof is much the same as in the small data case, with the pointwise lower 
and upper pointwise bounds coming from the functionals C[4>] and li[<j>], together with 
the norm bound H[4>], and an application of Proposition 14. 21 

We may prove global existence in H 1 as in the small data case by first showing 
that £ is conserved for H 1 solutions by taking a limit of H 2 solutions, 



Proposition 5.6 Let (j>o satisfy 



- 1 



< oo 



< oo 



(5.8) 
(5.9) 



and let </>(z, i; </>o) — 1 G C 1 ([0,T) : H 1 (M)) be the solution with data <f>Q. Then 
£[</>(•, t; 4>o)] is conserved on the interval [0, T). 

Proof: This proof is very similar to the small data case. Let e be the conserved lower 
bound, given by C[4>q] and let R = 1Z[4>q\ be the conserved upper bound. Requiring 
our strongly convergent sequence {<po, n } to satisfy 

||0o,7i - M\ H 2 < oo 



1 



< oo 



00, n 
ft[0O,n] < R 
£[0O,n] > e 

we proceed as in the small data case almost verbatim. 
Corollary 5.7 Theorem 15. 51 holds for H 1 . 



□ 



Corollary 5.8 For data satisfying and i)5.^|) . global solutions exist in C 1 ([0, oo) : 
H 1 (M.)) even when m = provided 4 < n < 6. 



We have thus established global well-posedness for data that may not be close to 
the uniform state in the H 1 norm, provided we satisfy I|5.1qI - I5.1c|I . 



6. Discussion 



Our definition of solutions of and our proofs of existence, both locally and globally 
in time, depend upon ensuring cf) > 0. It is this condition, and the use of Sobolev type 
inequalities to adhere to it, that limit our global existence proofs to the (n, m) pairs 
in the set defined by I|5.1al - I5.1c|) . However, numerical experiments indicate that we 
should expect a global lower bound on to exist at least for the values of [2, 5] x [0, 1]. 
On the other hand, solutions have been proposed in |141 1151 1131 which are zero on 
sets of measure zero. This leaves us with three questions: 
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Figure 4. The magma and generalized BBM equations, 16.11 and 16.21 

respectively, evolved from the data, 16.31 . initially bounded away from zero. While 
Magma remains bounded away from zero, BBM not only reaches zero but also 
becomes negative on a set of positive measure. 



• Why do solutions that are initially uniformly bounded below away from zero 
remain above zero? 

• How might a global existence proof be achieved for other choices of exponents n 
and m? 

• What is to be made of these functions that go to zero on sets of measure zero? 

We offer partial answers to these questions and suggest directions for future work in 
this final section. 

6.1. Importance of Nonlinear Dispersion for Positivity 

The appearance of a degenerate nonlincarity in the dispersive term appears to be key 
to the positivity of the solution. We saw this in Section [SJ where the appearance of 
terms oc ^ 2 ~™ - ™ 1 i n the invariants prevented <f> from going to zero. These terms only 
appear in the conserved integrals because of their presence in the underlying equation. 
There are physical reasons to suspect the benefit of a variable bulk viscosity in this 
model, consistent with (|5.1fe|l : by increasing m above zero, we gain global existence for 
a continuum of values of n, for which we could not otherwise succeed. But while having 
nonzero m may permit global existence for certain n, it is by no means required, as 
noted in Corollary 15. 81 

Rather than continue to look at the invariants, let us examine the evolution itself. 
Comparing 

<h + (^)z - (0 4 <M* = o (6.1) 

the magma equation with n = 4 and m = 0, shown to be globally well posed for 
arbitrary data in Sectional with a generalized BBM (gBBM) equation, 

u t + {u 4 ) x - u xxt = (6.2) 

we see a distinct and subtle difference in the dynamics; gBBM can become negative. 

With appropriate initial conditions, an assymetrical global minimum, 

/ ( x _ 75)2 \ 

no = O = 1 - .9exp (- j ^ _ — _ ^ j (6.3) 

(|6.2|) will cross the axis, while (|6.1|) will prop itself up above zero, as pictured Figure 
16.11 Away from the minimum in the figures, the solutions are quite similar. We 
suspect that it is this nonlinear dispersion that makes self-regularizing. 

6.2. A Comment on Invariant Bounds 

Better understanding the nonlinear dispersive term will likely be necessary to expand 
the existence result, as we believe the energy estimate approach, on its own, has 
been exhausted. Indeed, our estimates, such as (|5.3|l and l|5.4f) appear quite crude 
when evaluated numerically and compared to the maximum and minimum values a 
particular system will actually attain. 
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Figure 5. The evolution of data 16,41 for 16.11 . Note that it fails to become 
significantly larger pointwise from its initial maximum. 



Figure 6. The evolution of data 16.51 for 16.11 . Note that it fails to become 
smaller than its initial minimum of .5. 



If we start with initial data of the form 

M*) = 1 + (1-5 - 1) exp- ^" 1 ^ 2 (6.4) 

a Gaussian floating on a background of 1 with a peak at 1.5, (|6.1(l evolves into a 
solitary wave train, as in Figure^] The leading solitary wave has an amplitude on the 
same order as the data, about 1.6. If we were to use l|5.3[) to predict the maximum 
porosity, we would compute that U[<fio] ~ 10.6, an order of magnitude higher. 

Similarly, l|5.4|l predicts that the porosity will never go below a value many orders 
of magnitude smaller than the minimum value the profile eventually evolves into. For 
data for the form 

o (z) = l-(l-.5)cxp- (Z 2 ;^ )2 (6.5) 

corresponding to a negative Gaussian on a background of 1, with a minimum amplitude 
of .5, it never goes below .5. For the data in FigureEJ C[(j>o] ~ 8.1e — 05, not to mention 
that U[(j>o] ~ 26, again in excess of the prediction. 

6.3. Solutions that Reach Zero 

A question that remains is, can we start with smooth data bounded away from zero 
and reach zero? While no solution which starts above zero and subsequently reaches 
zero is known, several functions which are zero on sets of measure zero have been 
constructed. So called Compressive Solitary waves, with <f> < 1 and decay to 1 as 
£ — ► oc, are presented in |14U15j . Some of these are smooth and have a finite number 
of even roots. These are certainly classical solutions of (|l.l|l . although the conserved 
integrals of Section [JJ are not finite for them. 

One function in [TJ, along with some of those in [321 EH §° to zer0 at cus P s i 
behaving like 

4>(z,t) = 4(z - ct) <x y/\c\ 

near ( = 0. Such a function will not have a square integrable derivative, and 
consequently, as pointed out in for the case n — 3, m — 0, the conserved integrals 
will be also diverge. Thus these are not solutions in the sense of Definition 12.21 

However, they are solutions of the PDE in a weaker sense. We define a weak 
solution to (| 1 - 1 11 as a function <f>{z, t) such that 

-BtiKz, t)<t>(z, t) - d z ^{z, t)<t>{z, t) n + d z ^(z, t)4>(z, t) n d z ) dzdt = 0(6.6) 

for all ip(z, t) G Q°(R x M+), and such that both 
4>{z,t) 

'd t (j){z,ty 



(z,t) n d z 



( z ,ty 
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Figure 7. Figure (a) shows the phase portait for the traveling wave ODE. 16.81 . 
of 11.11 with c = 4. Figures (b) and (c) are weak solutions constructing using 
parts of the phase portrait. 



are in Lj oc (M.) in the z coordinate. Note that the solutions constructed in Section [5] 
will satisfy these conditions. 

Consider the case n = 3 and m = 0, with c > 3 and define F as follows 

{Curve A of Figure 7 (a) for Q < 
for C = (6.7) 

Curbe B of Figure 7 (a) for ( > 

Such a function is plotted in Figure 7 (b). For ( ^ 0, F solves the second order ODE 
obtained by making the traveling wave ansatz for and integrating up once with 
the boundary condition that F — > 1 at ±oo. 

-cF + F 3 + cF 3 F" = 1 - c (6.8) 

F is also continuous at C = 0. Expanding about ( = 0, 

F(C)~ V2(l-l/c) 1/4 VK\ 

and 

lim cF(C) 3 F"(C) = l-c 

If we define (cF 3 F")(0) = 1 - c, then (|6.7|) satisfies (|6.8|) as written, but not 
p" l - c 1 1 



cF 3 c F 2 
Since F satisifies (16 .8|) everywhere, 

0(z, t) 3 0(z, t) 2t = -c.F(z - ct) 3 F"(z -ct) = l-c- cF(z - ct) + F(z 

reducing the left-hand side of 1)6. 6fl to 



oo yoo 







(—d t ip(z, t) — cd z tp(z, t)) F(z — ct)dzdt 



This is just the weak form of the transport PDE 

d t u + cd z u = (6.9) 

which is solved by any function of the form u(z,t) — u(z — ct). Hence it is a weak 
solution in the the sense of (|6.6I) . It is bounded, continuous and in L i 1 oc (R). 
Unfortunately, there is a uniqueness problem. Let 

' Curve A of Figure E31» for C < 

for C = 

G(C) = < Curve C of Figure KHIla) for < C < z* (6.10) 

for C = z* 

^ Curve B of Figure E3» for ( > z* 

This function is pictured in Figure 7 (c). Except at and z*, G satisifes (|6.8(l pointwise. 
In addition, we claim that, as before, at the two cusps, 

lim cG(0 3 F"(0 = 1 - c 
|C|->o,**,c#>,»* 
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so that G in fact satisfies the ODE pointwise, everywhere. Using the same 
manipulation as above, we get 

POO />OG 

/ / (-d t t/}(z,t)-cd z ip(z,t))G(z-ct)dzdt 

JO J-oo 

Once again, we have reached the weak form of 1)6.9(1 and G is a weak solution to this, 
hence to the magma equation. Moreover, we could add any number of bumps of finite 
length of type C between the two curves A and B which go to 1 at ±00. 

A selection criterion is needed to distinguish between these functions in the 
same way an entropy condition is used to distinguish between physical and non- 
physical solutions to systems of conservation laws. Although even with such a 
selection criterion, we recall the physical argument from |12j that there will be stress 
singularities at a cusps, suggesting none of these are physically realizable. 
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